38.2 SCOPE vs STOP

Fit LME models for Prolactin, HPL and GH2.
Fit GEE models for Insulin and HOMA-IR.

• Outcomes: Prolactin, hPL, GH2, Insulin, HOMA-IR, Glucose (all outcomes log-transformed)
  
• Independent variables: Study, GA at sampling
  
• Model fitted: LME or GEE model
        log(hormones) ~ Study + GA at sampling + Study:GA at sampling
        LME only: Random intercepts for Study ID and Plate
        GEE only: Independence correlation structure

38.2.1 Global tests

Type III ANOVA on main and interaction effects, and LM model summaries.

PRL (ng/ml)

	Chisq	Df	Pr(>Chisq)
(Intercept)	2.473411	1	0.1	ns
Study	1.757843	1	0.2	ns
GA_samp	8.736573	1	0.003	**
Study:GA_samp	1.550230	1	0.2	ns

hPL (μg/ml)

	Chisq	Df	Pr(>Chisq)
(Intercept)	11.118300	1	0.0009	***
Study	5.570424	1	0.02	*
GA_samp	33.105635	1	<0.0001	***
Study:GA_samp	7.350297	1	0.007	**

GH2 (ng/ml)

	Chisq	Df	Pr(>Chisq)
(Intercept)	1.122593	1	0.3	ns
Study	5.234170	1	0.02	*
GA_samp	8.992223	1	0.003	**
Study:GA_samp	5.263896	1	0.02	*

Insulin (uU/ml)

	model term	df1	df2	F.ratio	Chisq	p.value
1	Study	1	Inf	0.395	0.395	0.5	ns
3	GA_samp	1	Inf	1.567	1.567	0.2	ns
2	Study:GA_samp	1	Inf	0.731	0.731	0.4	ns

HOMA-IR

	model term	df1	df2	F.ratio	Chisq	p.value
1	Study	1	Inf	0.298	0.298	0.6	ns
3	GA_samp	1	Inf	1.230	1.230	0.3	ns
2	Study:GA_samp	1	Inf	0.765	0.765	0.4	ns

Glucose (mmol/L)

	model term	df1	df2	F.ratio	Chisq	p.value
1	Study	1	Inf	0.006	0.006	0.9	ns
3	GA_samp	1	Inf	0.003	0.003	1.0	ns
2	Study:GA_samp	1	Inf	0.086	0.086	0.8	ns

38.2.2 Model summary

Ratio of Geometric means (95% CI)

PRL (ng/ml)

Var	Effect	P
(Intercept)	1.64 (-0.41, 3.67)	0.12
StudySTOP	1.46 (-0.69, 3.61)	0.19
GA_samp	0.20 (0.07, 0.33)	0.004
StudySTOP:GA_samp	-0.09 (-0.24, 0.05)	0.21

hPL (μg/ml)

Var	Effect	P
(Intercept)	1.82 (0.67, 2.96)	0.001
StudySTOP	-1.38 (-2.58, -0.17)	0.02
GA_samp	0.21 (0.13, 0.28)	<0.0001
StudySTOP:GA_samp	0.11 (0.03, 0.19)	0.008

GH2 (ng/ml)

Var	Effect	P
(Intercept)	-0.69 (-2.09, 0.78)	0.30
StudySTOP	-1.61 (-3.17, -0.09)	0.03
GA_samp	0.13 (0.03, 0.22)	0.006
StudySTOP:GA_samp	0.11 (0.01, 0.21)	0.03

Insulin (uU/ml)

Var	Effect	P
(Intercept)	1.37 (-2.34, 5.08)	0.47
StudySTOP	1.23 (-2.61, 5.07)	0.53
GA_samp	0.14 (-0.10, 0.38)	0.26
StudySTOP:GA_samp	-0.11 (-0.37, 0.14)	0.39

HOMA-IR

Var	Effect	P
(Intercept)	0.11 (-3.61, 3.83)	0.95
StudySTOP	1.08 (-2.79, 4.95)	0.58
GA_samp	0.13 (-0.11, 0.37)	0.29
StudySTOP:GA_samp	-0.12 (-0.38, 0.14)	0.38

Glucose (mmol/L)

Var	Effect	P
(Intercept)	1.71 (0.58, 2.85)	0.003
StudySTOP	-0.04 (-1.19, 1.10)	0.94
GA_samp	0.00 (-0.07, 0.08)	0.90
StudySTOP:GA_samp	-0.01 (-0.09, 0.07)	0.77

38.2.3 Diagnostics

Model diagnositcs.

Diagnositcs for PRL (ng/ml)

Diagnositcs for hPL (μg/ml)

Diagnositcs for GH2 (ng/ml)

Diagnositcs for Insulin (uU/ml)

Diagnositcs for HOMA-IR

Diagnositcs for Glucose (mmol/L)

38.2.4 Estimated means

Estimated marginal means (or adjusted means) adjusted for GA at sampling.

PRL (ng/ml)

Study	GA_samp	response	SE	df	lower.CL	upper.CL
SCOPE	10	37.84059	14.97325	183.71665	17.33445	82.60489
STOP	10	64.25777	10.79816	25.61462	45.47803	90.79244
SCOPE	16	125.40320	24.78532	46.73448	84.25636	186.64422
STOP	16	122.20695	26.05851	61.88348	79.79445	187.16261

hPL (μg/ml)

Study	GA_samp	response	SE	df	lower.CL	upper.CL
SCOPE	10	48.12118	9.400344	116.3177	32.68197	70.85399
STOP	10	35.67218	1.710800	43.9918	32.38566	39.29222
SCOPE	16	165.09165	11.552861	112.4283	143.71802	189.64396
STOP	16	234.19079	20.961034	179.3138	196.27551	279.43032

GH2 (ng/ml)

Study	GA_samp	response	SE	df	lower.CL	upper.CL
SCOPE	10	1.764828	0.4232374	27.87312	1.0797268	2.884636
STOP	10	1.052072	0.0640941	50.34783	0.9309196	1.188992
SCOPE	16	3.744028	0.3480694	126.03402	3.1148611	4.500281
STOP	16	4.304273	0.5007060	185.72738	3.4216229	5.414615

Insulin (uU/ml)

Study	GA_samp	response	SE	df	asymp.LCL	asymp.UCL
SCOPE	10	15.42851	10.443696	Inf	4.09391	58.14466
STOP	10	17.37172	1.623285	Inf	14.46449	20.86328
SCOPE	16	35.06535	5.621541	Inf	25.61042	48.01086
STOP	16	20.27501	4.212605	Inf	13.49286	30.46618

HOMA-IR

Study	GA_samp	response	SE	df	asymp.LCL	asymp.UCL
SCOPE	10	4.151566	2.7843932	Inf	1.115126	15.456106
STOP	10	3.826518	0.3878844	Inf	3.137036	4.667538
SCOPE	16	9.138816	1.7090499	Inf	6.334418	13.184787
STOP	16	4.200316	0.9587402	Inf	2.685292	6.570106

Glucose (mmol/L)

Study	GA_samp	response	SE	df	asymp.LCL	asymp.UCL
SCOPE	10	5.833361	1.0212224	Inf	4.139067	8.221199
STOP	10	4.947693	0.0651983	Inf	4.821543	5.077144
SCOPE	16	6.006065	0.4569785	Inf	5.173987	6.971957
STOP	16	4.740833	0.1711949	Inf	4.416896	5.088528

38.2.5 Comparisons

Hormone	Contrast	GA_samp	Ratio of Geometric means (95% CI)	P-value
PRL (ng/ml)	STOP / SCOPE	10	1.70 (0.81, 3.58)	0.2
PRL (ng/ml)	STOP / SCOPE	16	0.97 (0.67, 1.42)	0.9
hPL (μg/ml)	STOP / SCOPE	10	0.74 (0.50, 1.10)	0.1
hPL (μg/ml)	STOP / SCOPE	16	1.42 (1.15, 1.75)	0.001
GH2 (ng/ml)	STOP / SCOPE	10	0.60 (0.36, 0.99)	0.05
GH2 (ng/ml)	STOP / SCOPE	16	1.15 (0.87, 1.51)	0.3
Insulin (uU/ml)	STOP / SCOPE	10	1.13 (0.30, 4.30)	0.9
Insulin (uU/ml)	STOP / SCOPE	16	0.58 (0.35, 0.97)	0.04
HOMA-IR	STOP / SCOPE	10	0.92 (0.24, 3.48)	0.9
HOMA-IR	STOP / SCOPE	16	0.46 (0.26, 0.82)	0.008
Glucose (mmol/L)	STOP / SCOPE	10	0.85 (0.60, 1.20)	0.3
Glucose (mmol/L)	STOP / SCOPE	16	0.79 (0.67, 0.93)	0.005

38.2.6 Barplots by Study

Barplots in this section may be misleading due to involvement of Study by GA interaction.

38.2.7 Plots for GA